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Abstract 

We present a practical S'-matrix to potential inversion procedure for coupled- 
channel scattering. The inversion technique developed is applied to non- 
diagonal Sill ^^"^ SY>ui one projectiles, yielding a tensor interaction Tr and is 
also applicable to spin-i plus spin-i scattering. The method is a generaliza- 
tion of the iterative-perturbative, IP, method. It is tested and evaluated and 
we investigate the degree of uniqueness of the potential, particularly for cases 
where there is insufficient information to define the potential uniquely. We 
examine the potentials which result when the S'-matrix is generated from a 
Tp interaction. We also develop the generalisation, using established proce- 
dures, of IP S-matrix-to-potential inversion to direct observable-to-potential 
inversion. This 'direct inversion' procedure is demonstrated to be an efficient 
method for finding a multi-component potential including a Tr interaction 
fitting multi-energy o", iTn, T20, T21 and T22 data for the scattering of spin-1 
nuclei from spin-zero target. It is applicable to other channel spin 1 cases. 

PACS numbers: 21.30.-x, 13.75. Cs, 25.10.-hs 
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I. INTRODUCTION 



Various methods for carrying out ^-matrix to potential inversion are now available, 
see for example [^0, but, until recently, it has been possible only for cases with channel- 
spin zero or 1/2. However, there have been many accurate experiments involving spin-one 
polarised particles and these provide a powerful motivation to develop an efficient technique 
for inversion in cases with higher channel spin, i.e. coupled-channel scattering. With such a 
technique, one can exploit the very large volume of polarisation data which has accumulated. 
This includes vector and tensor analysing powers and polarisation transfer observables. 

Over the years, we have developed a practical and widely generalisable procedure, the 
iterative perturbative, IP, method and we recently [|] demonstrated an extension to 
spin-1 projectiles for the first time. Ref. |^ demonstrates the method in one specific appli- 
cation, but does not give details or a derivation of the method. In this paper we present 
details of our extension of the IP S'-matrix to potential inversion method to the coupled- 
channel case of spin-1 projectiles and present further evaluation of it. We also test and 
evaluate an important extension of the IP method, single step data-to-potential inversion 
for coupled- channel scattering. 

In Ref. P] we applied single step inversion to analyse real data. We note here that there 
are also many ways in which IP S ^ V inversion can contribute to understanding nucleus- 
nucleus interactions. Perhaps the most obvious application is the inversion of S'-matrix 
elements found by phase shift analysis of experimental data. However, there are also many 
important applications which involve the inversion of theoretical S'-matrix elements, i.e. 
elastic channel S derived from coupled channel calculations, Glauber model and resonating 
group model calculations. The potential found in this way contains information concerning 
the contribution of tensor components to dynamic polarization and exchange processes to 
inter-nuclear potentials. For the particular case of spin-polarised deuteron and ®Li scattering, 
obvious applications include the study of the infiuence of reaction channels and distortion 
effects on the projectile-nucleus potential, especially in its non-central components. There 
are a number of longstanding puzzles relating to spin-polarised deuteron scattering, including 
the anomalously small real part of the tensor interaction, which can be studied using these 
methods. 

Coupled-channel inversion represents a significant development in inversion techniques 
since a non-diagonal potential is derived from a non-diagonal S-matrix. Such a potential 
couples channels of the same conserved quantum numbers but different values of orbital 
angular momentum. Spin-1 inversion is therefore the first example of coupled channel in- 
version in which a non-diagonal S'-matrix is made to yield a non-diagonal potential. Apart 
from the general derivation, the present paper is framed for a very specific two-channel case: 
deuterons scattering from a spin zero nucleus with inversion determining a tensor interac- 
tion involving the non-diagonal operator Tr as given below. In spite of the rather specific 
nature of the present application, we believe this work opens the way to a fully general class 
of coupled-channel inversion situations involving the determination of a coupling potential 
from a non-diagonal S'-matrix. 

This paper presents in detail only those aspects of the IP inversion formalism which 
are connected with the specific coupled channel generalisation to spin-1 scattering. In other 
respects it calls upon previous publications 1^-0 in which the general aspects of the IP 
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inversion procedure are described. Because of the specific application to spin-1 projectiles, 
we establish our notation by beginning in Section II with a brief review of basic aspects of 
spin-1 scattering. 

An important feature of the IP S ^ V inversion procedure is the natural way in which 
it can be convoluted with (data) —>■ S fitting to give an overall (data) —>■ V algorithm. 
This provides a new and efficient data analysis tool which in many cases obviates the need 
for independent (data) S inversion. Important advantages follow when fitting data for 
many energies since the underlying potential model guarantees that the energy dependence 
of the S-matrix will be smooth without the need to postulate parameterized forms for S{E). 
Indeed, we have shown [y] that (data) V inversion can provide a powerful alternative 
method for phase shift fitting of multi-energy scattering data for light nuclei. 

Ref. contained a restricted analysis of low energy d + ^He data. A future paper will 
present a much more exhaustive analysis of the very large collection of data for this system. 
At a later stage we hope to present an analysis of ^Li+'^He data, including tensor analysing 
powers. 



II. THE SCATTERING OF SPIN-1 NUCLEI FROM SPIN-0 TARGETS 
A. Formalism for spin-1 scattering 

In order to establish our notation, we outline the standard formalism for the elastic 
scattering of spin-1 projectiles from a spin-0 target in the presence of tensor forces. The key 
feature introduced by the tensor interaction of Tr type (see below for a classification of tensor 
forces) is that it couples channels of different orbital angular momentum /. Specifically, for 
particular values of the conserved quantities J, the total angular momentum, and vr, the 
parity, whenever vr = (— l)"'"'"^, then two values of orbital angular momentum, I = J — 1 and 
J + 1, are coupled by Tr. 

For total angular momentum J and orbital angular momentum /' the radial wavefunction 
ipin satisfies the coupled equations. 



^, + k^-'pJ\v\i'J) 



+ 1) 



where /i is the reduced mass of the system and {V J\V\l" J), a function of r, is the matrix 
element of the inter-nuclear interaction V integrated over all angular and internal degrees of 
freedom. The second subscript, I, on ip identifies the incoming orbital angular momentum. 
This is determined by imposing on the solution of the coupled equations, the following 
asymptotic boundary conditions : 



^f,,ik,r)-.6inIi'ikr)-S-^iOiikr) 



(2) 



Here, Ii{r) and Oi{r) are the incoming and outgoing asymptotic Coulomb radial wavefunc- 
tions, often written Hi{r)* and Hi{r) respectively as in Satchler |]TD|], namely: 



Iiikr) = Gi{kr) - iFi{kr)- Ii{kr) = Gi{kr) - iFi{kr) 
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where Fi and Gi are regular and irregular Coulomb wavefunctions respectively. Note that 
the boundary conditions given in Eq. ^ differ by a factor from those adopted by Satchler. 
Where there is no ambiguity, we suppress the J superscript. When vr = (— 1)"^, V is diagonal 
and Eq. |l] is uncoupled. 

In general, S will not be unitary, but will be subject to the unitarity limits: l^np + 
1 5*121^ < 1 and 15*22^ + |5'2ip < 1, where, of course Su = 821- These limits present no 
particular problem for S V inversion where S can be assumed to satisfy them, but they 
can represent a significant problem in the case of data to potential inversion, see Section [V^ . 



B. The Tr interaction and its effect. 

Non-diagonal matrix elements {l'J\V\l"J) occur in Eqn. |I]for elastic scattering of spin-1 
projectiles with certain types of tensor force. The possible forms of local tensor interaction 
have been classified by Satchler [|T^,0 who defined Tr, Tl and Tp interactions. The Tl 
interaction is believed |jT^ to be very small, at least below 50 MeV/u, and is in any case 
diagonal in /. The Tp interaction could well be substantial JT^I but appears to be hard to 
distinguish phenomenologically from Tr. The gradient operators within Tp make calculations 
harder, and the present inversion method does not apply to it. 



1. The Tr operator 



In this work we assume that the inter-nucleus potential V contains a tensor force com- 
ponent of Tr form [|r^ : 



Tj,Vnir) = {{s-Ty-2/3)Vn{r). 



(3) 



We quote the matrix elements of the interaction Tr for future reference. The diagonal 
matrix elements of Tr are: 





1 = J-l 


/ = J 


/ = J+ 1 


< JI\Tk 


Jl > 


1 J-l 

3 2J+1 


1 

3 


1 J +2 
3 2J+1 



and the non-diagonal matrix elements are: 



< J J - 1|Tr|JJ + 1 >=< J J + 1|Tr|JJ - 1 >= 



2J+ 1 



(4) 



2. The radial form of the Tr interaction 

The derivation of VR(r) from the folding model has been discussed at length long ago 
by Keaton and his colleagues [p!6H l8|| and also by Raynal Within the folding model, 
the deuteron Tr interaction arises directly from the D-state component. The overall gen- 
eral success of folding models for central and spin-orbit interactions suggests that folding 
model calculations of VR{r) should give at least approximately the correct radial form and 
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overall magnitude, but this has not been borne out in the case of Tr according to extensive 
phenomenological studies, e.g. The overall conclusion is that the real part of VR(r) 



predicted by the folding model is much too strong for heavy nuclei, of the right order of 
magnitude for light target nuclei, and actually about three times too strong for a "^He target. [| 
These facts, together with a large literature discussing breakup and reaction channel con- 
tributions, suggest that we have no generally applicable reliable knowledge of VR(r). There 
is reason to doubt even the general arguments, based on folding models, that it should be 
small in the interior of heavier nuclei, away from the nuclear density gradients in the surface. 
Such gradients define the angle between the projectile spin s and the vectorial position r of 
the projectile with respect to the nuclear centre, see Eq. ^. 

III. S -^V INVERSION FOR SPIN-1 PROJECTILES ON SPIN-0 TARGETS 

A. General background of IP inversion 

The IP method has been successful for Sij — > V{r) + I ■ sVso{r) inversion for spin half 
projectiles, and we now present its generalisation to spin-1 inversion. The only restriction is 
to a Tr tensor interaction. Certain features of the IP method, to our knowledge not shared 
by other inversion methods, will be of particular importance in the particular systems to 
which we shall apply spin-1 inversion. These include the ability to find an explicitly energy 
dependent potential from phase shifts for a range of energies, the ability to handle a range 
of energies simultaneously and to include Majorana terms for all potential components. For 
many applications the important property is that mentioned in the introduction, i.e. that 
IP inversion lends itself to direct observable to potential inversion. This not only avoids 
the need for independently determined phase shifts (or S'-matrix), but actually provides an 
advantageous method of determining such phase shifts. For a full description of IP inversion 
as applied in the spin- 1/2 case see Refs. 1^-0]. The formalism presented in Ref. 0, whereby 
energy dependent potentials are obtained from multi-energy datasets, can be used with spin- 
1 inversion as described here, although energy dependence is not actually exploited in the 
test cases. A brief general account of IP inversion is given in the next section. 

B. IP inversion for the coupled channel case; application to spin-1 

Our notation must reflect the fact that the outcome of inversion will be a potential 
with many components. We therefore label each component with an index p which identifies 
central, spin-orbit or tensor terms, each real or imaginary. The number of components 
doubles when the potential is parity dependent. (Parity dependence is particularly important 
for light nuclei at lower energies.) 



^Ref ||8| exploits the inversion formalism presented here to give an alternative analysis of d + ^He 
scattering, a theme elaborated in later papers. 
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The IP method commences with a 'starting reference potential', SRP, and proceeds by 
iteratively correcting each component p of the potentiah 

V(P) -> r + J2 (^n'^^n^ (r) (5) 
n 

where a^^^ are coefficients to be determined and vl^\r) are the functions comprising the 'in- 
version basis', (which, if required, can be chosen differently for different p). The amplitudes 
a^^-* are determined at each iteration from linear equations, based on an SVD algorithm, 
which successively reduce the 'phase shift distance' a defined by: 



cr 



'-T.\si-si\'. (6) 



For each partial wave k, is the 'target' S-matrix and 5*^ is for the potential at the 
current iteration. Here the label is a single index which identifies the partial wave angular 
momentum / as well as the energy Ei when multi energy sets of Si{Ei) are simultaneously 
inverted. It also includes non-diagonal elements of Sfi, in the spin-1 case described later. 
The linear equations are based on the (usually) very linear response P,^, AS", of the 



complex S'-matrix to small changes AV^ in the potential. The expression for this is well 
known in the uncoupled case and is very simple: 

^Si = T2T mr)f^V{r)d.r. (7) 



h'kJo 

In Eq.|^, the S'-matrix Si is written in terms of the asymptotic form of the regular radial wave 
function as 'ipiir) Ii{r) — SiOiij-) where Ii and Oi are incoming and outgoing Coulomb 
wave functions as before. When inverting Si{Ek) over a series of energies Ek-, the energy label 
Ek is implicit in these equations with index k subsumed with orbital angular momentum 
I to give an overall channel label. In the case of spin-1/2, the j label is also subsumed in 
the same way ||^. Linear algebraic equations for local variations of a^^-* follow from the 
minimisation of cr^, [^^. 

We now present the generalized linear response relationship which applies to the non- 
diagonal S-matrix for spin-1 elastic scattering. The derivation is given in Section [111 B 1| . 
For any given set of conserved quantum numbers, certain channels will be coupled by the 
nucleus-nucleus interaction and we use labels k. A, /x, v for these channels. Thus the matrix 
element of the nucleus-nucleus interaction V between the wavefunctions for channels n and 
A, corresponding to integrating over all coordinates but r, will be written Vi^\{r). The 
increment AS^a in the non-diagonal S- matrix which is due to a small perturbation AVka('"), 
is 

^S.x = ^J2 r i^^..ir)AV,,ir)i;,xdr (8) 
h k Jo 

where ipt,,^ is the z/th channel (first index) component of that coupled channel solution for 
the unperturbed non-diagonal potential for which there is in-going fiux in channel k (second 
index) only. The normalisation is ip^^^ — > — Su^Oi^ where /; and Oi are incoming 

and outgoing Coulomb wavefunctions for orbital angular momentum /; there is no complex 
conjugation in the integral. Starting from Eq.H, spin-one inversion becomes a straightforward 
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generalisation of the procedure outlined above and described in Refs. [|5|-0. The method is 
implemented in the code IMAGO [26| where the linearity relations have been exhaustively 
tested by the gradient method. 

A convenient feature of the IP method is that one can judge from the behaviour of o"^ 
as the iteration proceeds whether a satisfactory inversion has been achieved. A low value of 
cr^ obviously guarantees that a potential closely reproducing the input Si has been found. 
Because the IP method is implemented interactively, there is an opportunity to examine the 
potential for oscillatory features. These might well be spurious and result from over-fitting 
noisy data. In such a case, one can reduce the basis dimensionality or raise the SVD limit 
and this generally allows one to achieve a smooth potential, often with only a small increase 
in a^. One must bear in mind that genuine oscillatory features, corresponding to non- locality 
in an L-independent local potential or to L-dependence of the underlying potential, can be 
necessary to achieve a precise representation of Sij or S/q. 

We stress here that, although the context of our discussion is the determination of 
a tensor interaction from the non-diagonal S'-matrix elements describing the scattering of 
spin-1 projectiles, the range of application is much more general. 



1. Derivation of non- diagonal perturbation expression. 

We now outline the derivation of Eq.|]. The derivation can be applied to the general 
coupled channel inversion from non-diagonal S'-matrix to non-diagonal potential. Our start- 
ing point is Eq.|l| which we shall write with a simplified notation for two channels. Until the 
last step in the argument, we shall assume we are using units in which ti^ /2m = 1. 

The radial wavefunctions in channel i with incoming waves in channel A are written as 
ipix and have asymptotic behaviour at r ^ oo: 

i^ix Sixh - SxiOi (9) 

where for simplicity we write Oi for the outgoing Coulomb wavefunction with orbital angular 
momentum ^ appropriate to channel and similarly for the ingoing wavefunction /j. For 
brevity we omit labels for conserved quantum numbers J and vr. 

Absorbing the centrifugal interaction within the potential, we can write the coupled 
equations for the radial wavefunctions appropriate to incoming waves in channel A as: 

€a = T.(y^, - ES,j)^jX, t = l,2. (10) 
j 

In the case considered here, Vij with i ^ j arises entirely from the tensor interaction, 
the matrix Vij being symmetric. Now, denoting by if) the wavefunction arising from a 
(symmetric) perturbation in the potential Vij —>■ Vij + AVij, we can write: 

^^X = E(^U- + ^y^J - E5,,)i,,x. ^ = 1, 2 

i 

Multiplying Eq.[l^ by Tpi^j, and Eq. 
one from the first one, we get: 
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(11) 



11 by ipia, summing over i and then subtracting the second 



E 3Z (^a^.M - ^'^X^^,) = - E A V^^.^.A, (12) 

as all terms including Vij in the right hand side vanish due to symmetry of Vij. 

Integrating Eq. ([l^) from r = to the asymptotic region and using the usual Wronskian 
relationship W[l£,Oe] = —2ik, we get: 

_ roo _ 

2ik{S^X - 5a^) = E / i^^Ay^3i'3>^^^■ (13) 

Since the coupled channel equations with a symmetrical potential matrix give a symmetrical 
^-matrix = S^\, the left hand side of Eg . (|T3|) is equal 

2ik{S^x - S^x) = 2ikAS^x (14) 

Hence we find, reinstating the 2m/ fi^ factor on the right hand side: 

iui - 
^Sx, = ^EX V^^AAl^,,%^dr. (15) 

This expression is valid for any symmetrical finite perturbations AVij decreasing at r — > cxo 
sufficiently rapidly. It is easily extended to any system of coupled-channel equations. For 
small AVij we make the Born approximation assumption that ip ^ if) and get: 

A^Am = t27:E/ i^^xmJi^J,dr. (16) 
n k^ Jo 

The expression (|T6|) is the basis for the coupled-channel inversion method. The success of 
the inversion method in leading to a converged solution confirms the wide applicability of 
this equation in each step of our iteration process. 



IV. TESTING COUPLE-CHANNEL S INVERSION FOR SPIN-1 

PROJECTILES 



We carried out two contrasting tests of S* — inversion as described in Sections [IV B 
and [IV C| below. First let us define the potentials and the inversion basis. 



A. Specification of the interaction potential and basis used 

In one respect our notation is non-standard. We write down the complete potential for 
spin-1 projectiles scattering from a spin-zero target. It is 

Ken(r) + il^cen(r) + Koui(r) + 21 ■ s(V;o + iW^so) + (^R + iWR)Tn (17) 

where KouK'^) is the usual hard-sphere Coulomb potential. Note that our spin-orbit poten- 
tials Vso and Wso, are defined in such a way that they will be half the magnitude of those 
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defined according to the usual convention for spin-1 projectile.0 For tlie test cases we 



present, tlie spin-orbit potential is defined as in Eq. [17. 

For notational simplicity, Eq. |1^ has not been written to reflect parity dependence. 
There are two alternative methods of representing parity dependence. The code IMAGO 
can apply either of these to each of the components in Eq. |1^ except Koui- The first 
representation defines Wigner and Majorana components for each term, say 

K = K,w + (-l)Vx,M (18) 

where I is the partial wave angular momentum. With this form the inversion procedure 
can be made to determine V^,w and V^^u for any or all V^. An alternative approach is to 
determine independent positive or negative parity components for V^. In many cases, the 



Wigner- Majorana representation is most natural and has been shown [|ri|] to be preferable 
where the odd parity term may otherwise be ill-determined. However, sometimes the odd- 
even representation is more appropriate, for example where a particular has completely 
different shapes and magnitudes for the different parities, as we believe can be the case 
for Vr. The code IMAGO offers the freedom to represent the parity dependence of each 
component in either way. 

The IP method is not tied to any particular set of functions for the inversion basis and 
each component of the potential can be represented by a different basis. It is an important 
feature of the IP method, as implemented in IMAGO, that a range of different functions 
is available, and those which we have applied are specified in [^. Zeroth order Bessel 



functions and harmonic oscillator functions are both linearly independent sets which have 
proven useful where bases of large dimensionality are necessary to describe a potential over a 
wide radial range down to r = 0. For cases involving light nuclei, particularly for inversion of 
small S'-matrix datasets, a small basis comprising a series of Gaussian functions is preferable. 
A Gaussian basis covering just the nuclear surface region is also useful for heavy ion cases 
where there is no information available to determine the potential in the nuclear interior. 
It is important that a basis should not be chosen which would describe the potential over a 
radial range, or to a degree of detail, which is not warranted by the information contained 
in the set {Si} or by the nature of the physical situation. In practice, much smaller bases 
are often necessary in order to eliminate spurious oscillatory features from the potentials. 
The operation of the SVD algorithm, with adjustable SVD limit, stabilises the inversion and 
can, where appropriate, reduce the effective dimensionality of the inversion basis. 

B. Single-energy inversion 

In Ref. we presented a test of S ^ V for deuterons scattering from the hght nucleus 
'^He in which there is very little absorption. Here we present a test for a much heavier nucleus, 
where there is substantial absorption, and demonstrate that a potential, very accurate almost 
to the nuclear centre, can be obtained by inversion. The test case studied was for a ^^Ni 



^In our papers relating to spin 1/2 projectiles, tlie usual convention has been used. 
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target and deuterons at a laboratory energy of 56 MeV. The parameters found by Hatanaka 
et al fitting angular distributions and the analyzing powers Ay and Ayy were used. A 
notable feature of the potential was that the imaginary Tr term was quite large, but the real 
Tr term was very small (a common but unexplained feature of deuteron optical potentials.) 
The spin-orbit component was real. The potential was parity independent as expected for 
this combination of target, projectile and energy. The energy and other characteristics of the 
reaction are such that there are 'many' active partial waves. 'Many' here means sufficient, 
even with for a single energy, to yield a precise reproduction of the potential. 

The test was carried out as follows: one of the authors applied the optical model pa- 
rameters of Hatanaka et al to the standard spin-1 scattering code DDTP reading out 
the S'-matrix onto a file. A second author, knowing only the target and the energy, then 
applied IMAGO to find the potential from these S'-matrix elements. The inversion was car- 
ried out with a starting potential which contained only real and imaginary central potential 
components. These were guessed from general systematics without specific knowledge of 
Hatanaka's potential. Since with IMAGO there is complete freedom to choose the starting 
potential and inversion basis, it is worthwhile to test the inversion method starting with no 
more information about the potential than might be available in a 'for real' case. 

When a converged solution was found, the potentials obtained were compared with the 
known potentials with results shown in Figure 1. The solid lines represent the 'target' po- 
tential, i.e. the potential from which the S'-matrix was calculated using DDTP. The short 
dashes represent the potential found by inversion; it can be seen to reproduce the target po- 
tential very closely except very near to the nuclear centre in the case of the imaginary tensor 
component. We have not shown the real tensor component which was only a few percent of 
the imaginary component nearly everywhere. This small component was reproduced only 
qualitatively, as expected, since the absolute errors for the real and imaginary Tr compo- 
nents were similar in magnitude and comparable to the real Tr potential itself. In Figure 
1, the dashed line represents the starting potential, zero for the spin-orbit and tensor com- 
ponents. The S'-matrix elements for the target and inverted potentials are indistinguishable 
on a graph, corresponding to values of a of roughly 10~^. 

This test shows that the inversion procedure has the capability of revealing quite fine 
details of the potential as would be required for the kind of applications, discussed in the in- 
troduction requiring the inversion of single energy S derived from theory. Such studies might 
establish, for example, the contribution of specific exchange terms or reaction couplings to 
the inter-nucleus potential. 

C. Multi-energy inversion at very low energy 

At low energies and for light target nuclei, very few partial waves are involved so that 
there will in general be insufficient information contained in the S'-matrix elements for a 
single energy to yield a detailed and precise potential. The situation is even worse in cases 
where parity dependence must be assumed since this halves the information available for 
potential components of each parity. The problem can be ameliorated if S is available for 
more than one energy. If S is available over a narrow range of energies, then the algorithm 
can be made to yield to an energy independent potential; this is what we have earlier called 
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'mixed case' inversion (see Refs. and first of Ref. 0]) and, in effect, the information 
from the energy derivative of S is exploited. In many cases, Sij or Sfq are provided over 
a wide range of energies. In this case one should ideally consider the potential to be en- 
ergy dependent and determine the energy dependence itself. This can be done within the 
framework of the parameterisations presented above. 

An example of where the sets of 5*^^ or S/q are too small to define the potential very 
closely is the S ^ V inversion situation embedded in the analysis of low energy, experi- 
mentally determined, multi-energy observables for d -|- alpha scattering. A first report was 
presented in Ref. ||^. The test we now describe is directly relevant and asks the following 
question: what properties of the potential can reliably be determined from very small sets 
of ^? 

The test was for deuterons scattering from ^He with Sfn calculated from a known poten- 
tial at 11 energies: 8, 8.5 . . . 12.5, 13 MeV. The known potential was energy independent but 
parity dependent and was taken to be real. (The imaginary parts of empirical potentials are 
known to be small for d + '^He at these energies.) The following terms were included: central 
Wigner, central Majorana, spin-orbit Wigner and separate even parity and odd parity Tr 
tensor potentials (the odd/even choice for Tr refiects what we believe P to be the case for 
the actual d + ^He tensor force.) The central and spin-orbit terms are like those found in 



Ref. |H, and the very large even parity tensor term is based on that of Dubovichenko [29 



see also 111 



The inversion was effectively 'mixed case' in the sense just described. The starting 
potential was zero in all components except for the Wigner real central and Wigner real 
spin-orbit terms. In keeping with the nature of this test, the very small inversion basis 
of Ref. [§] was used. This has two Gaussians only for each component except the central 
components for which there were three. The centres and widths of the Gaussians were not 
varied during the inversion. 

The 'target' (known) and inverted potentials are shown in Figure 2, together with the 
starting potential required by the IP method. The starting potential is the dot dashed line, 
non-zero for two components only, and corresponding to a = 10.552 where a, defined in 
Section [111 B| , is summed over the 11 energies. The inverted potential is shown as the dotted 
line, and the 'target' potential, from which 5*/^ was calculated, is the full line. We see that 
the qualitative features are reproduced although less well for the small components and near 
the nuclear centre. The value of a for the potential shown in the dotted line was 0.135 which 
is reasonable for a low energy multi-energy case. For 10.5 MeV, this corresponds to S/q for 
the target and inverted potentials being indistinguishable on a graph apart from one single 
term: the phase angle of the non-diagonal part of Sfq for higher partial waves for which, 
in any case, the magnitude {Sfql, I ^ I', is very small. The tensor potential, having very 
different odd and even parity components, is as well reproduced as could be expected with 
the very small basis. Note that the starting potential for the inversion has zero tensor terms. 
From the matrix elements of Tr given in Section |II B| , we see that I = partial waves are 
ineffective and hence we cannot expect to reproduce the tensor real term ^ at r = 0. 

In summary: we found that the qualitative properties of the potential were reliably 
reproduced, particularly for the larger components. Thus, reliable statements about the 
general features of d -|- *^He potentials can be made, but nothing can be asserted concerning 
non-central interactions for r < 0.5 fm. 
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V. INVERTING S/,^ CALCULATED WITH A Tp TENSOR INTERACTION. 

The inversion technique which we have described is hmited to a tensor force of the Tr 
type. Since there exist processes which are expected to lead to Tp forces, the possibihty 
must be faced that data analysed using the data-to-^ extension of the inversion method, 
which is described in Section will indeed involve a Tp tensor interaction. It is therefore 
relevant to ask, in the context of S" — inversion: can we invert Sf^ calculated with a Tp 
tensor interaction with a potential which has only a Tr tensor interaction? If so, to what 
extent does inversion yield valid central and spin-orbit components? 

There is further interest in knowing how well the general effects of a Tp interaction can 
be represented by a Tp potential. The properties and even existence of a Tp interaction 
have not yet been convincingly linked to experiment since the consequences of the two 
kinds of interaction are difficult to distinguish phenomenologically. This was discussed by 
Goddard who compared Sf,i and the observables calculated from a Tp interaction with 
the corresponding quantities calculated from a particular Tp interaction devised in such a 
way that, according to semi-classical arguments, it would be very similar in effect. 

We study these questions by exploiting the equivalent pairs of tensor potentials intro- 
duced by Goddard. We first inverted S/q for 30 MeV deuterons scattered from ^^Fe with a 
Tp potential and then inverted S/n derived from the potential containing that Tp interaction 



which is 'equivalent' in Goddard's sense. The two potentials are given in Table 1 of Ref. [30 



The first part of the test showed that inversion of for a known Tp still works very well 



at about half the energy of the test described in Section [IV B| . The results were very similar: 
the Tp potential, which in this case is of a volume Woods-Saxon form with depth 5 MeV, is 
accurately reproduced even at the nuclear centre. The solid and (scarcely distinguishable) 
dashed lines in Figure 3 respectively represent Goddard's original potential and that found 
by inversion. The Sfq for the inverted potentials, including the non-diagonal terms, are 
indistinguishable on a graph from those for the original potentials. 

The dotted lines in Figure 3, show the inversion for Goddard's Tp case. The non-tensor 
components are qualitatively well reproduced, the derived potentials having the appearance 
of the target potentials but with superimposed oscillations. This waviness is relatively more 
significant for the small components, the real central potential being reproduced to within 
a few percent for all r. The Tp interaction found by inversion is now surface peaked in 
form but of average depth comparable to that of the Woods-Saxon (which however had a 
local momentum dependence, see |^|). The diagonal S/q for target and inverted potentials 



are graphically indistinguishable, as are aigS/q for / 7^ /' for low values of J. However 
the non-diagonal ^-matrix was not well reproduced for J > 7, for which partial waves the 
non-diagonal \Sf,i\ is very small. The value of a was much higher than for the Tp case, i.e. 
0.0294 compared with 0.00589. 

The results presented graphically in Figure 3 can be quantified in terms of the volume 
integrals and rms radii for the central and spin-orbit components of the inverted potentials. 
For the Tp case, all of these quantities were reproduced to a few parts in a thousand with 
the (small) volume integral of the spin-orbit term being least accurate: the error was 0.7 
%. The errors for the non-tensor components found when inverting Goddard's Tp potential 
were a few percent, the volume integral of the spin-orbit term again being least accurate 
with an error of 3.8 %. 
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Goddard also performed an identical comparison for the case of 13.0 MeV deuterons 
scattering from ^^Ti, and we repeated the test just described for this case. There is interest 
in doing this since the inversion algorithm applied to S/q for a single energy is expected to fail 
at lower energies for reasons explained in Section |IV C| . However, we find that the results for 



13 MeV deuterons on ^^Ti are essentially the same as for 30 MeV deuterons on ^^Fe for both 
Tr and Tp interactions. The form of the Tr potential representing the actual Tp component 
was essentially the same as that shown for 30 MeV in the bottom panel of Figure 3 and this 
similarity applies also to the deviations of the non-tensor terms. It therefore appears that 
we have found general properties of the Tr potential representing an actual Tp potential. 

As a result of these tests, and noting that Tp interactions are not predicted to be 
particularly large, we conclude: 

1. The existence of processes of the kind which give rise to a Tp component will not 
prevent this inversion procedure, which includes only Tr tensor interactions, from 
fitting S/,i and is unlikely to greatly falsify inversions of this kind, particularly with 
regard to the non-tensor components. IP spin-1 inversion as described here is thus 
not fatally undermined by the possible existence of Tp interactions. The effort needed 
to develop spin-1 inversion including Tp interactions require greater motivation than 
exists at present. 

2. As Goddard suggested, almost all the effects of such a potential can be well represented 
by a Tr tensor interaction, although its relationship to the form of the Tp interaction 
is, as might be expected, more complicated than can be deduced from simple semi- 
classical arguments ||3^. The phenomenological problem of establishing Tp interactions 



is still considerable. 



VI. DATA TO POTENTIAL INVERSION FOR SPIN-1 PROJECTILES 

In what follows, we first briefly indicate how S —>■ V inversion is extended to (data) —>■ V 
inversion for the uncoupled case, then indicate how this is extended to include coupling, as 
is required for spin-1 scattering. 



A. Data to potential inversion for uncoupled situation 

For clarity we suppress spin-related subscripts and begin by recasting Equation |^, using 
Equation |, as @J^J^: 

f^Qj 1777 poo 

' ' (^z(r))2i;(P)(r)dr. (19) 



We now introduce a conventional function: 



x'-Y.{-^] (20) 

k=l \ ) n k=l\ ^^kn J 
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where a^^ and Pj^ are the input experimental values of cross sections and analyzing powers 
of type n (cr, iTn, etc.) respectively. When fitting data for many energies at once, the index 
k indicates the angle and also the energy. Data re-normalising factors can be introduced as 



an additional contribution to Equation ^ 



We must now expand in terms of the q;^^\ To do this we first linearize the calculated 
cross sections and analyzing powers, by expanding ak (and Pkn) about some current point 
{a(f)(z)} (see Ref. §): 




Aal^\ (21) 



which applies at each iterative step i = 0, 1, 2,. . . and the correction (to be determined) for 
the n-th amplitude ). Equivalent relations are applied for the P's. 

Linear equations result from demanding that is locally stationary with respect to 
variations in the potential coefficients a^^-*, i.e. the derivatives of with respect to the 
potential components a^^^ must vanish. Solving these linear equations is straightforward for 
any reasonable number of them and yields corrected values al^\i). We then iterate the whole 



procedure, with wave-functions ipi in Equation |T9| calculated using the corrected potentials 
from Equation ^, until convergence is reached. This algorithm almost always converges 
very rapidly, in general diverging only when highly inconsistent or erroneous data have been 
used or when the iterative process involves a very unsuitable starting point. Multi-energy 
(data) — > V inversion is thus reduced to the solution of simultaneous equations in a series 
of iterative steps. 



B. Generalisation to spin-1 



Spin-1 (data) V inversion is a natural generalisation of the above formalism with Si 
replaced by S/q and Eq. |19| replaced by the analogous form derived from Eq.|[ It is shown 
in Ref. ||^ that the system does indeed converge to a potential which fits the observables. 



C. Evaluation of ambiguities for spin-1 (data) V inversion 



The tests of (data) V inversion must refiect the way it will be applied; this is rather 
different than for S ^ V inversion. With the latter, one often has quite precise S calculated 
from a theory, and one then seeks quite precise and subtle properties of V, often relating 
to modifications of the theory. Inversion from measured observables is different because the 
data are generally far from complete and will contain statistical and, possibly, systematic 
errors. For this reason, we must be less ambitious concerning the details of the potential to 
be extracted. The test therefore ask the following question: for a situation with few active 
partial waves, how well-determined can we expect the potential to be? 

As in 5* — > \^ inversion, one must never attempt to establish details of the potential for 
which the input data carries no information. We must therefore apply the smallest possible 
inversion bases and accept approximate solutions. The penalty for excessive inversion basis 
dimensionality is the occurrence of spurious oscillatory features. In effect, at low energies 
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where the data is incomplete and featureless (reflecting the small number of partial waves), 
the goal of (data) V inversion is to flnd the smoothest potential compatible with the 
data. IP inversion affords a level of control in this respect that is not possible with other 
inversion procedures. 

The test we describe is for low energy d + '^He scattering. The results will be useful for 
interpreting previous flts to experimental data for this system. The following observables, 
a, iTii, T20, T21 and T22, were calculated at laboratory energies of 8, 9, 10, 11, 12 and 13 
MeV using the same purely real, energy independent potential used in Section |IV Q Apart 
perhaps from the extremely strong 'Dubovichenko-type' tensor interaction, very strongly 
peaked at r = 0, the general features of this potential are, we believe, similar to those 
of potentials which flt actual experimental data. This energy range is somewhat above 
the broad 2+ resonances and the region of strong mixing between the l"*" channels. The 
observables were evaluated for the six energies over a range of 20*^ to 170° CM, at intervals 
of one degree, and Gaussian noise was added as follows. For a, 1% errors were imposed. For 
iTii, the errors were 2% of the maximum magnitude and for the three tensor observables, 
5% of the maximum magnitude. 

We then applied (data) —>■ V inversion to this multi-energy dataset, seeking a single 
energy independent potential. Following Section [IV (J| and Ref. the inversion bases for the 
Wigner and Majorana real central components consisted of three Gaussian functions. For the 
other components there were just two Gaussians. As in Section |1V Q , the starting potential 
was zero in all but the Wigner real central and Wigner real spin-orbit components. The 
results are shown in Figure 4 where we compare the known ('target') potential (solid lines), 
the chosen starting potential of the iterative method (dash-dotted lines, two components 
only), and two inverted potentials, shown as dashed and dotted lines. The dashed hues 
show the potential found after a flrst sequence of iterations and correspond to x^/^ = 15.473 
where F, the number of degrees of freedom, was ~ 4500. This number arises since we seek 
simultaneous flts to flve observables at six energies and 150 angles. The effective number 
of parameters is ~twelve. At this stage the reproduction of the larger components of the 
potential is fair, but the tensor terms are poor, with the even parity real tensor term being 
still almost zero. The corresponding flt to the model data is indicated by the set of dashed 
lines in Figure 5. The flt is of a quality which would be widely regarded as quite good when 
fltting experimental data, with only T22, and perhaps T20 around 120 degrees, fltted poorly. 
The quality of flt to T21 is remarkable in view of the very poor reproduction of the tensor 
interaction. 

A subsequent further set of iterations led to an almost perfect flt with / F = 1.2155. 
Figure 4 shows that the potential, dotted, flts all parts of the potential except at quite small 
radii. In particular, the even parity real tensor is perfectly fltted for r > 1 but not fltted 
at all for r < 1. This is in accord with arguments given in Section [IV Q . As expected from 
the values of jF, the flts to the 10 MeV dataset, shown as dotted lines in Figure 5, are 
essentially perfect, being scarcely visible over the angular range of the artiflcial data. The 
same potential simultaneously flts the observables for the other flve energies comparably well. 
We conclude that we could not expect to establish the various components of the potential 
to a higher degree of accuracy than shown in Figure 4 by fltting available experimental 
data. It is very salutary to see, in Figure 4, the profound change in the nature of the tensor 
interaction which follows the improvement of the flt revealed in Figure 5, comparing dashed 
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and dotted lines. The intermediate inversion, dashed hnes in Figure 5, represents a fit of a 
quahty which is often deemed acceptible when fitting experimental data. We note without 
further comment that the desirability of pursuing the best possible phenomenological fits is 
sometimes called into question. 

It should be noted that the computing time required on a modern workstation to carry 
out the direct inversion of the data is very modest, and certainly much less than required 
to carry out a model independent optical model search, particularly one involving odd and 
even parity Tr components and about 4500 degrees of freedom. 

In Ref. we discussed the application of direct inversion of data as a method for phase 
shift analysis. It is therefore of interest to see the quality of fit to 5*/^ which corresponds 
to the two fits shown in Figure 5. The top three panels of Figure 6 show the phase shifts 
corresponding to the l = J — 1, l = J, l = J + l diagonal components of S, and the 
bottom panel presents half the argument of the non-diagonal S. The solid lines show the 
known potential, the dashed line is for the /N = 15.473 fit and the solid line is for 
the x^/N = 1.2155 fit. For two of the panels, the solid and dash-dot lines are nearly 
indistinguishable but they are clearly distinguishable in the other two, suggesting that there 
are limits to phase shift determination even when over some 4000 data are fitted with 
X'^/N = 1.2155. 

We conclude that direct inversion is a practical, reliable and efficient means of de- 
termining a local potential which represents large, multi-energy datasets including tensor 
observables. The example presented here indicates the extent to which the results obtained 
by this method are meaningful at low energies where few partial waves are involved. 

VII. SUMMARY AND CONCLUSIONS; SURVEY OF POSSIBLE 

APPLICATIONS 

We have presented details of an inversion procedure which can be applied both to 
spin-1 projectiles scattering from a spin-0 target nucleus and to spin-| plus spin-| particle 
scattering. The non-diagonal Sfq yield a non-diagonal potential containing a tensor term. 
To our knowledge, this is the first time this has been achieved, and opens up the possibility 
of wide range of other inversion scenarios ranging from other channel spin-1 cases (such as 
p + scattering) to the inversion of S'-sub-matrices of higher dimensionality. There are 
many other capabilities inherent in the underlying IP method: these include the possibility of 
inverting Sfq for several energies leading directly to an energy dependent potential, including 
bound state energies within the input data, and the ability to handle cases where parity 
dependence must be allowed for. 

In this paper we have presented tests of IP S* ^ V spin-1 inversion and evaluated its 
performance in 'difficult' cases. We showed that when there are sufficient active partial 
waves, the procedure yields very accurate potentials even quite near the nuclear centre. 
Where, on the other hand, there are few partial waves available to define each potential 
component, as is typical with light nuclei at low energies and where the potential is parity 
dependent, it is still possible to extract the qualitative features of a potential. 

We also addressed the fact that the method is at present limited to Tr tensor interactions 
although it is quite probable that processes leading to Tp interactions are active. We showed 
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that Sf,i arising from Tp interactions can be fitted quite well with a Tr tensor interaction 
and that, moreover, this does not lead to serious errors in the non-tensor components of the 
potential. 

The IP inversion algorithm also forms the basis of a very efficient alternative way to find 
a multi-component local potential which fits elastic scattering data, particularly for multi- 
energy datasets. This is the direct (observable) V inversion procedure in which the IP 
S ^ V inversion is embedded. This 'direct inversion' can be applied to spin-1 projectiles. 
We examined the ambiguity problems which arise in a 'difficult' (i.e. few partial waves, 
parity dependence) test case which is relevant to the evaluation of an analysis of low energy 
d + "^He scattering, the subject of a recent and an extended future publication. Known 
potentials can be very well re-fitted, but it is clear that the non-central terms cannot be well 
established at the nuclear centre. In the course of performing this inversion test, it became 
apparent that fits of widely accepted quality lead to tensor potentials which have nothing 
in common with those determined by pursuing 'perfect fits'. 

Finally, we remark that the method we have demonstrated here is certainly not limited 
in usefulness to deuteron scattering. It would certainly be worthwhile applying it the elastic 
scattering data for halo nuclei when these are of sufficiently substantial information content. 
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FIGURES 



FIG. 1. Potential for deuterons scattering from Ni at 56 MeV. Prom top, the real central, 
imaginary central, real spin-orbit and imaginary tensor terms. Solid lines represent the target (i.e. 
input) potential, dots the potential found by inversion and the dash-dot the starting potential of 
the inversion. The starting potential was zero for the real spin-orbit and the tensor terms. There 
was no imaginary spin-orbit potential; see text regarding the real tensor term. 

PIG. 2. Por deuterons scattering from ^He for energies from 8-13 MeV, the starting (dash-dot), 
inverted (dotted), and known, 'target', (solid) potentials. Prom the top, the potential components 
are: Wigner real central, Wigner real spin-orbit, even parity real Tr, Majorana real central and 
odd-parity real Tr. 

FIG. 3. For deuterons scattering from ^^Fe at a laboratory energy of 30 MeV, comparing 
components (from top, the real central, imaginary central, real spin-orbit and real Tr) of the 
potential of Goddard (solid line) having a Tr term, the potential (dashes) found by inverting Sf^ 
calculated from the solid line potential and the potential (dotted) found by inverting Sfn calculated 
from Goddard's 'equivalent' Tp-containing potential. The dashed and solid lines are only clearly 
distinguishable for the smallest, i.e. spin-orbit, component. 

FIG. 4. For deuterons scattering from '^He for energies from 8-13 MeV, potentials found by 
direct inversion of artificial data incorporating Gaussian noise as described in the text. The starting 
potential, two components only, is shown as a dash-dotted line, the true ('target') potential is solid 
and the first stage {x^/N = 15.473) and second stage {x^/N = 1.2155) inversion potentials are 
shown dashed and dotted respectively. 

PIG. 5. Por deuterons scattering from ^He for energies from 8-13 MeV, the fits to the 10 

MeV part of the six energy artificially generated datasets. The final fit with /N = 1.2155 is the 
dotted line (visible in a few places only) and the intermediate fit with 'x^ /N = 15.473 is shown 
dashed. 

PIG. 6. Por deuterons scattering from ^He for energies from 8-13 MeV, fits to the diagonal (top 
three panels) and non-diagonal ^-matrix. The solid curves correspond to the known potential, the 
S'-matrix for the /N = 15.473 fit is dashed and for the /N = 1.2155 fit is dotted. As described 
in the text, the four panels are essentially the (real) phase shifts. 
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